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Abstract 

We review a new iterative procedure to solve the low-lying states of the 
Schroedinger equation, done in collaboration with Richard Friedberg. For the 
groundstate energy, the n^^ order iterative energy is bounded by a finite limit, 
independent of n; thereby it avoids some of the inherent difficulties faced by 
the usual perturbative series expansions. For a fairly large class of problems, 
this new procedure can be proved to give convergent iterative solutions. These 
convergent solutions include the long standing difficult problem of a quartic 
potential with either symmetric or asymmetric minima. 
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1. Importance of Low- lying States 

It is a common belief that most of natural phenomena should be described by 
solutions of the Schroedinger equation. Take the nonrelativistic Schroedinger 
equation of a system of charged particles with Coulomb interactions. In prin- 
ciple, this single N-dimensional second order linear partial differential equation 
contains in its solutions all information about many disciplines: atomic and 
molecular physics (except for small relativistic corrections), condensed matter 
physics, chemistry and biology. Furthermore, for most of these applications 
we need only knowledge about its low-lying states. For these low-lying states, 
actually we do have very good information about their wave functions when 
the Coulomb interaction is strong; i.e., when the wave amplitude is large. How 
come that we still have great difficulties in solving these equations? Quite 
often, the challenge lies in how to handle the large number of configurations 
when the interaction is weak. 

In quantum chromodynamics (QCD) and quantum electrodynamics (QED) 
the tasks are even more difficult since we do not have any exact solution. 
Besides the very promising lattice and other numerical calculations, we rely 
mostly on perturbative expansions. Yet, such expansion often leads to a diver- 
gent series with zero radius of convergence, as exemplified by the perturbative 
calculation of the gyromagnetic ratio g of the muon[l-3]: 

= -- + 0.765857376(-)2 

Z Z TT TT 

+24.05050898(-)^ + 131.0(-)^ 

TT TT 

+677(40)(-)^ + --- (1.1) 

TT 

where a is the fine-structure constant. The coefficients are all positive, and 
each successive one becomes larger and larger. There are good reasons [4] to 
believe that the vacuum state of QED would be unstable against pair creations 
if a could be analytically continued to the negative region. Thus, the difficulty 
of perturbative series seems to be again closely related to our inability to solve 
the low-lying states of the corresponding Schroedinger equation. In QCD, 
because of the existence of instanton configurations, similar serious problems 
also exist for its vacuum and other low-lying states. 
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This situation may be illustrated by the following one dimensional problem 
with a quartic potential of degenerate minima. The Schroedinger equation can 
be written as 



where T = and 



{T + V - E)iIj = (1.2) 



V{x) = \g\x' - a'f. (1.3) 



An alternative form of the same problem may be obtained by setting q = 
^/2ga{a — x), so that the Hamiltonian becomes 

H = T + V = 2gan (1.4) 

with 

in which H, q and e = l/^/8ga^ are all dimensionless. The perturbative 
expression of the groundstate energy E is[5-14] 

1 9 1 89 1 f)( ^ ) 

or, in terms of the dimensionless coupling 

E-2aal\-e^--le*-^p -■■■], (1.7) 

which has a similar characteristics as the {g — 2) expansion ( 11 .11) . 

In this one-dimensional case, from ( 11 .51) . one sees that for negative, H 
has no groundstate; this explains why the power series ( II .71) is divergent. For 

= l/8ga^ positive, it is also possible to trace the origin of the difficulty. 
The corresponding perturbative series of the groundstate wave function il){x) 
may be written as 
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Since the potential V{x) has two degenerate minima at x = ±a, there are two 
different expansions depending on the normalization condition: 



i/jia) = 1, or i/ji-a) = 1. (1.9) 

In either case, the perturbative series is divergent for any g. However, for g 
large, the wave function i(j{x) can be approximately described by the first two 
terms of these expansions. For x > 0, using the normalization i/j{a) = 1, we 
find 

a + x ' 

and for X < with tlj{—a) = 1, 

a — X 

Unlike ^(x), neither expression has a zero derivative at a: = 0. Furthermore, 
each carries a spurious pole-term, atx = —a for (11.101) and a: = a for (11.111 ). In 
addition, if one stays with the perturbative series (11.81) . each successive higher 
order expansion accentuates further these spurious pole terms. This suggests 
that we must not follow the high order perturbative expressions, especially in 
the region x near 0, when i/j{x) is exponentially small. 

In all these problems, it is not difficult to construct a trial function such 
that 

(f) K, ifj when iIj is large, (1-12) 

and with the correct symmetry properties of il) (e.g., 0(x) = 0(— x) for the 
quartic potential problem (ll.2l) -( |T751) ). The question is how to device a con- 
vergent procedure that can lead from the trial function to the correct ijj. In 
the following we will discuss such a method. 
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2. The New Method 



As examples of the types of Schroedinger equations that we are interested, 
consider the bosonic component of QED or QCD in the axial gauge, or a system 
of nonrelativistic particles. In either of these cases the kinetic energy T is a 
quadratic function of momentum operators. Through a linear transformation 
of the relevant coordinate variables, the Hamiltonian H can be written as 

H = T + V{q) (2.1) 

where 

Q. = {Qi, Q2, • • • , Qn) (2.2) 

and 

1 N 1 

T=-^ZdVdq^ = --V', (2.3) 

^ i=l ^ 

with = oo in the case of a field theory. Let ip^q) be the groundstate 
determined by 

{H-E)i/; = 0. (2.4) 

To derive ip{q), we proceed as follows [15-19]: 

1. Construct a good trial function 0(g), so that 

c/>{q) « ij{q) (2.5) 

in regions when ^(g) is expected to be large. An example of how to construct 
such a trial function is given in Appendix A. 

2. Construct U{q) — Eq by differentiating (f){q): 

U{q) -E,^ (^VV)/0. (2.6) 



Define 



H^ = -\^^ + U{q). (2.7) 
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We have 

{Ho - Eo)cl){q) = (2.8) 

where the constant Eq may be defined by requiring, say, the potential energy 
part U{q) of Hq to satisfy at g = oo 

U{oo) = 0. (2.9) 

Alternatively, in problems like the quartic potential ( 11 .31) . the constant Eq can 
be more conveniently determined by requiring 

minimum of U{q) = 0. (2.10) 

(Since our method deals only with the operator Hq — Eq, it is independent of 
the particular way to define £"0.) Introduce the differences h and S by 

HQ-H=U{q)-Viq) = h{q) (2.11) 

and 

Eq-E = E. (2.12) 
The original Schroedinger equation ( 12.41) can be written as 



{HQ-EQ)^{q) = {h-8)^{q). (2.13) 
For the discussions of the groundstate, we take both 

il){q) and 4){q) to be real and positive. (2-14) 



Multiplying ( 12.131 ) by 0(g) and ( 12.81) by ij{q), then taking their difference, 
we have 

- • (0V^ - V^V0) = {h- S)(t) il). (2.15) 



The integration of ( 12.151) over all space gives !{h — E)4>iljd^q = 0; therefore. 



S = J h (j) d^q / J (pip d^q. (2.16) 



It is convenient to introduce 

/(g) = V^(g)Mg); (2.17) 

(12.161 ) becomes then 

S = l h{q)ct>\q)f{q)d\ / / ct?{q)f{q)d''q. (2.18) 



3. We propose to solve the Schroedinger equation ( 12.131) by an iterative pro- 
cess, defined by 

{H, - E,)^,,{q) = {h{q) - Sn)^n-i{q) (2.19) 
with n = 1, 2, • • • and for n = 0, 

Mq) = Hq)- (2.20) 



Multiplying ( 12.191) by (f){q) and ( 12.81) by i/Jn{q), then taking their difference, we 
derive 

- • {(j)VllJn - ^nV0) = {h- 8n)(t> V^n-l- (2.21) 



Analogous to (12.161) . we have 



Sn = jh(j) A / / A- (2.22) 

Note that in each step of the iteration from n — 1 to n, for any particular 
solution i^n^q) that satisfies ( 12.191) . so is 



Mq) = ^l{q) ^ cnct^iq) (2.23) 

also a solution where c„ is an arbitrary constant. Since 0(g) is positive in 
accordance with ( 12.141) . we can choose so that 



il)n{q) is positive everywhere. (2.24) 

Define 

fn{q) = Mq)/Hq) (2.25) 



and therefore ( 12.221) becomes 

Sn = Jh (j)^fn-id''q j J (f^fn^id^'q. (2.26) 
From ( 12.241) . we have 

fn{q) positive. (2.27) 

As we shall see, a particularly convenient choice of the constant c„ is to set 
the minimum of fnio) to be 1 and therefore 

fn{q) > 1. (2.28) 



From ( 12.261) and ( 12.271) . it follows that is bounded if h{q) is bounded. 



Therefore, as n ^ oo, limSn cannot be oo. This avoids the type of divergence 
difficulties that appears in the example (11.71) . It is therefore reasonable to 
expect that lim„^oo should converge to the correct S, at least when h 
is small. As we shall see, for a fairly large class of problems, including the 
quartic potential (11.31) . the convergence of £„, turns out to be independent 
of the magnitude of h{q), provided that it is finite and satisfies some general 
conditions. 

In this paper, we concentrate on the groundstate. Extensions to some low- 
lying states can be found in Ref.[17]. 
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3. An Electrostatic Analog 

In terms of the ratio fn = introduced in (12 .2 51) . the n^'^ order iterative 

equation ( 12.211) can be written as 

- • = an{q) (3.1) 



where 



and when n = 0, 



^n{q) = {h{q) - Sr,)4>\q)fn-M (3.2) 



k{q) = 1 (3.3) 

in accordance with (12.201) . Assuming that i/Jn-iiq) has already been solved, 
we can determine through (12.221) . Therefore, a^iq) is a known function. 
Consider a dielectric medium with a q-dependent dielectric constant, given 

by 

K{q) = c^\q). (3.4) 

Interpret cr„(g) as an external electrostatic charge distribution, i/„ the elec- 
trostatic potential, —\Vfn the electrostatic field and 

Dn = -\^Vfn (3.5) 
the corresponding displacement vector field. Thus (13.11) becomes 

V ■ Dn = cr„, (3.6) 

the Maxwell equation for this electrostatic analog problem. 

At infinity, 0(oo) = 0. In accordance with (13.41) - (13.51) . we also have 
Dn{oo) = 0. Hence the integration of (13.61) leads to the total external elec- 
trostatic charge to be also zero; i.e., 

/ an{q) d^'q = (3.7) 
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which is the same result given by ( 12 .261) for the determination of £n- Because 
the dielectric constant ti{q) in this analog problem is zero at = oo, the di- 
electric media becomes a perfect dia-electric at oo. Thus, the equation of zero 
total charge, given by ( 13 .71) . may serve as a much simplified model of charge 
confinement, analogous to color confinement in quantum chromodynamics. 
We note that (13.11 ) can also be derived from a minimal principle by defining 

I{Uq)) = / {\^{yfnf - CInfn}d''q. (3.8) 

Because of ( 13 .71) . the functional I{fn{q)) is invariant under 

fniq) fn{q) + constant. (3.9) 

Since the quadratic part of I{fn{q)) is the integral of the positive definite 
^/^(V/„)^, the curvature of I{fn{q)) in the functional space fn{q) is always 
positive. Hence, I{fn{q)) has a minimum, and that minimum determines a 
unique electrostatic field — |V/„, as we shall see. To establish the uniqueness, 
let us assume two different V/„, both satisfy ( 13.11 ). with the same k = (j)'^ 
and the same cr„; their difference would then satisfy ( 13.11) with a zero external 
charge distribution. For cr„ = 0, the minimum of I{fn{q)) is clearly zero with 
the corresponding V/„ = 0. To derive fn{q) from V/„, there remains an 
additive constant at each iteration, as already noted in ( 12.231) . As we shall 
discuss in the next section, this arbitrariness allows us the freedom to derive 
different types of convergent series. 

We note that in this electrostatic analog, the problem depends only on two 
input-functions: 0^(g) and the product 

hiq)(PHq). (3.10) 

The original potential V{q) no longer appears explicitly. Consider, e.g., a 
system of particles with electric charges. For configurations when the Coulomb 
energy V{q) is attractive and near its singularity, we have good knowledge 
about the wave amplitude ^^{q), which is large but regular. The same should 
be true for a good trial function 0(g). In accordance with ( 11.121) . we expect 

h{q) small when i/j is large. (3.11) 
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For configurations when i/j is small, so should be 0. Thus, a good trial function 
(j) would result in a small everywhere even though h may not. 

When = 1, we may denote the coordinate ( 12.21 ) by a single x. Corre- 
spondingly f l3.5lj and f l3.6lj become 

Dn = -l^f^ (3.12) 
2 ax 



and 

dPn 
dx 

Since Dn{oo) = 0, we have 



= an{x). (3.13) 



Dr,{x) = - J an{z)dz (3.14) 

a; 

and therefore 

oo oo 

Ux) = /,(oo) - 2 / r'(?/)c^2/ / <Jn{z)dz (3.15) 



with 



which satisfies 



an{x) = {h{x) - 8n)4>^{x)fn-i{x), (3.16) 



oo 

j an{x)dx = 0. (3.17) 



Therefore, 



D,(-oo) = D„(oo) = (3.18) 

and 

oo oo 

/„(oo) - /„(-oo) = 2 / (p-\y)dy J an{z)dz. (3.19) 

-oo y 

When N > 1, for radially symmetric problems we have 

h = h{r), (j)'^ = K{r) (3.20) 
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where 

r' = g? + gi + --- + 4- (3.21) 

Likewise i/j, i/jn and cr„ are all functions of r only. The radial component of 

Dn is 



(D„}r = -l<r)^. (3.22) 



Correspondingly, ( 13.61) becomes 

^|:(-"-'{A.),)=-n{r), (3.23) 

with 

an{r) = {h{r)-6n)(t>\r)fn-i{r). 

Hence, 

-j oo 

{Dn)r = -pvTT / z^'-'aniz) dz (3.24) 



r 



and on account of ( 13.221) 



00 



Mr) = /..(OD) - / / z"-'ajz) dz. (3.25) 

Assuming as r ^ oo, 

02(r) ^ a e-^"' (3.26) 
with a, b and / all positive constants. We assume 

/ > 2. (3.27) 
In accordance with ( 13.161) . as r ^ oo 

an{r) c e-^'' (3.28) 

where 

c = (/i(oo) - Sn) a /„_i(oo). (3.29) 
12 



By partial integrations, we see that as r ^ oo, 

fn{r) - /n(oc) - - _ ^y^, [1 + ^(^)J • (3.30) 

When the dimension = 1, the same expression applies, with r replaced by 

X. 

Because of (13.71 ). in the > 1 radially symmetric case we have 

00 

iV-1 



r^"V„(r) dr = 0; (3.31) 
therefore ( 13.251 ) is identical to 







/n(r) = /n(0) - / -T-T-IrT / ^^"V.(z) dz. (3.32) 



2 dy ^ 

When N = 1 and V{x) is an even function of x, ( 13.320 applies with the 
replacement of r by x; if V{x) is not an even function then we have, from 
dHHD and dUl), 

X V 

fn{x) = /n(-oo) - 2 y (j)'^{y)dy J an{z) dz. (3.33) 



-oo 



For an arbitrary > 1 dimensional radially non-symmetric problem, the 
general solution of the electrostatic analog problem can no longer be reduced 
to simple quadratures. Thus, numerical analysis may become important. In 
that case, the minimal principle using the functional I{f{q)) given by ( 13.81) 
may be of some practical use. 
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4. Hierarchy Theorem 

In this section we restrict our discussions to either (i) the A^^-dimensional radially 
symmetric case in which the functions V in ( 12 . 11) and h in ( 12.111 ) are 

V = V{r) and h = h{r) (4.1) 



where r'^ = qf + ql + ■ ■ ■ + q% as \r\ ( 13.211) . or (ii) the one-dimensional case 
( I3.12I) -( ITT7D with even functions 



V{x) = V{-x) and h{x) = h{-x). (4.2) 
In the latter case, we need only to consider the region 

x = r>0. (4.3) 

Thus, we need only to consider (i), since the A^-dimensional radially symmetric 
solution ( 13.251) reduces to the one-dimensional case ( 13.151) . with = 1 and r 



replaced by x. Furthermore, we assume h{r) to satisfy at all finite r > 0, 

h{r) > and h'{r) < 0; (4.4) 

at infinity, 

h{oo) = 0. (4.5) 

Throughout the paper, ' denotes ^. An example of such h{r) is given in 
Appendix A for the quartic potential problem (II. 3j) . 

In accordance with ( 13.251) . each n*^ iterative solution fn{r) carries an in- 
dependent additive constant. In the following, we differentiate two sets of 
sequences, labelled A and B, satisfying different boundary conditions: 

/„(oo) = 1 for all n, in Case (^4) (4.6) 

or 

/„(0) = 1 for all n, in Case (B). (4.7) 
Thus, by using ( 13.251) or ( 13.321) we have in Case (^4) 

fnir) = 1-2 r y-''+'r\y)dy T z^-V'(^)(/^(^) " S,,)^.,{z)dz, {4.8A) 

Jr Jy 
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and in Case (B) 



' z''-^4>\z){h{z) - 8n)fn-l{z)dz. (4.8S) 



In both cases, En is determined by the corresponding /n-i(r) through ( 12.221) 
i.e., 



in which [F] of any function F{r) is defined to be 

[F] = r r^-^(P^{r)F{r)dr. 

J 

It is convenient to introduce x(^) defined by 
Thus, ( 14.101) becomes 

[F] = I x\r)F{r)dr. 

Let Pn{f), and Vn{r) be related to cr„(r), «:(r) and (i?ri)r of (13.21) . 
and (p^ by 



(4.9) 
4.10) 
4.11) 
4.12) 



K,{r) = ^K{r) = x^(^) 



and 



P.(r)=r^-^(DO, = --/C(r)/;(r) 



4.13) 



4.14) 



4.15) 



Correspondingly, in Case (^4), (4.8A) can be written as 



fn{r) = 1 - 2 X-\y)dy P x\z){h{z) - Sn)fn-i{z)dz, {4.16A) 



and in Case (B), (4.8B) becomes 



fn{r) = 1-2/^ X 



X^{z){h{z) - Sn)fn~l{z)dz. 



(4.165) 
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Because of ( 13 .2 31) . we have 

V'^{T) = pn{T) (4.17) 

and therefore 

Vn{r) = - / pn{z)dz = - / {h{z) - £n)x\z)fn-i{z)dz. (4.18) 

Jj' Jr 



From (B7|l . 

I Pn{r)dr = (4.19) 

which leads to 

1^n{r) = r Pn{z)dz = fjhiz) - 8n)x\z)fn-l{z)dz. (4.20) 



These two expressions of Vn{r), (I4.18j ) and (I4.20j) . are valid for both cases 
(^4) and {B). Let be defined by 



h{r) — 8n = Q at r = 



(4.21) 



Since h'{r) < 0, ( 14.211) has one and only one solution, with h{r) — Sn negative 
for r > Tn and positive for r < r„. Thus, if 



fn-i{r) > 
for all r > 0, we have from i KW\ and i K7^ 

Vn{r) > 
and therefore, on account of ( 14. 151) . 

m < 0. 

In Case (^4), because of /n(oo) = 1, ( 14.241 ) leads to 

/n(0) > fn{r) > fn{00) = 1. 



(4.22) 



(4.23) 



(4.24) 



(4.25^) 



Since for n = 0, fo{r) = 1, ( 14.221 ) - (4.25A) are valid for n = 1; by induction 
these expressions also hold for all n; in Case (^4), their validity imposes no 
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restriction on the nnagnitude of h{r). In Case {B) we assume h{r) to be not 
too large, so that (4.16B) is consistent with 

/n(r) > for all r 

and therefore 

/,(0) = 1 > fnir) > /,(oo) > 0. (4.255) 

As we shall see, these two boundary conditions (A) and (B) produce sequences 
that have very different behavior. Yet, they also share a number of common 
properties. 

Hierarchy Theorem (A) With the boundary condition /n(oo) = 1, we have 
for all n 

Sn+l > Sn (4.26) 

and 

Thus, the sequences {£n\ and {/n(^)} are all monotonic, with 

8i<£2<Sz<--- (4.28) 

and 

l</i(r)</2(r)</3(r)<--- (4.29) 

at all finite r. 

{B) With the boundary condition /n(0) = 1, we have for all odd n = 2m+l 
an ascending sequence 

^1 < ^3 < ^5 < • • • , (4.30) 
but for all even n = 2m, a descending sequence 

e2>e4>S6>---; (4.31) 
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furthermore, between any even n = 2m and any odd n = 2/ + 1 

> S21+1. (4.32) 
Likewise, at any r, for any even n = 2m 



dr \ f2m{r) 
whereas for any odd n = 2/ + 1 

d (f2i+2{r] 



< 0, (4.33) 



dr \}2i+i{r)) 

The proof of the hierarchy theorem is given in Appendix B. The hierar- 
chy theorem was discussed in our earlier papers for one-dimensional prob- 
lems[18,19]. The extension to the A^-dimensional radially symmetric case is 
new. As we shall also discuss in Appendix B, the lowest eigenvalue E of the 
Hamiltonian T + F is the limit of the sequence {En} with 

En = Eo- 8n. (4.35) 

Thus, the boundary condition fn{oo) = 1 yields a sequence, in accordance 

with iKm . 

Ei> E2> E2,> ■■■> E, (4.36) 

with each member En an upper bound of E, similar to the usual variational 
method. 

On the other hand, with the boundary condition /,i(0) = 1, while the 
sequence of its odd members n = 21 + 1 yields a similar one, like ( 14.361) . with 

Ei> E2,> E^> ■■■> E, (4.37) 

its even members n = 2m satisfy 

E2 < E^< Eg < ■■ ■ < E. (4.38) 

It is unusual to have an iterative sequence of lower bounds of the eigenvalue 
E. Together, these sequences may be quite efficient to pinpoint the limiting 

E. 

In Appendix C, we discuss the comparison between the iterative solution 
and perturbative series for some simple examples. 
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5. Tribute 



This paper is based on my joint work with Richard Friedberg. It was a pleasure 
for nne to present the main body of the paper at the Number and Nature 
Symposium held on December 17-18, 2004 at the Rockefeller University, as a 
celebration for Mitchell J. Feigenbaum's 60th birthday. In the Chinese tradition, 
sixty is very special. It signifies the successful completion of one's first cycle, 
and the fresh beginning of a new second cycle. 

Because China has no indigenous religion, there is no origin of reference in 
the Chinese calender, nothing equivalent to BC or AD. All Chinese years are 
labelled cyclically, modulus sixty. 

In Chinese calender, each year is designated by two indices. The first index 
runs cyclically from one to ten, with each number represented by a different 
word shown in the first column of Chinese characters in Figure 1. We may 
translate the first character at the top of that column as alpha (for one) and 
the second character as beta (for two), then gamma (for three), • • • until kappa 
(for ten); after kappa, we have alpha, beta, • • • again. The second column of 
Chinese characters has twelve words, each for a different animal. In Figure 
1, the year 1944, when Mitchell was born is designated (Alpha-Monkey), the 
next year 1945 is (Beta-Chicken), • • •; ten years later, the first index for 1954 
becomes again alpha, but the second index is horse; two years later for 1956, 
the first (number) index is gamma, but the second (animal) index becomes 
monkey again. It takes | • 10 • 12 = 60 years to complete the cycle. This year 
2004 is like 1944, the Alpha-Monkey year and Mitchell is now in his new cycle, 
starting a new beginning. 

the most beloved story in China is the book "Monkey" by Wu Cheng-En 
written in the 16th century. It is about an immortal Monkey-King of Chaos. 
The translation was by Arthur Waley (Figure 2) in 1943, perhaps anticipating 
that 1944 was the beginning of Mitchell's Alpha-Monkey cycle. To celebrate 
2004, the beginning of Mitchell's newest cycle. Figure 3 is a drawing of mine, 
entitled " Monkey King the Thinker" . All of us are looking forward to this new 
cycle with our very best wishes. 
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Figure Captions 

1. Ciiinese calender is cyclic, modulus sixty. 

2. One of the nnost beloved novels in China 

3. " Monkey King the Thinker" by T. D. Lee (2004) (color online) 
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Chinese Calender 



(gamma) i\ (dog J 



(kappa ) ^ 



^ t 



(sheep) 
(fnonkey) 



Y f (montey) 
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Monkey 
King of Chaos 

Written by Wu Cheng-En 
(16th Century) 

Translated by Arthur Waley 
(1943) 



Figure 2 
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Monkey King The Thinker 



Figure 3 
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Appendix A Construction of Trial Functions 
A.l A New Formulation of Perturbative Expansion 

As mentioned in Section 1, in many problems of interest, perturbative ex- 
pansion lead to asymptotic series. Nevertheless, the first few terms of such an 
expansion could provide important insight to what a good trial function might 
be. For our purpose, a particularly convenient way is to follow the method 
developed in Refs.[15] and [16]. As we shall see, in this new method to each 
order of the perturbation, the wave function is always expressible in terms of a 
single line-integral in the N-dimensional coordinate space, which can be readily 
used for the construction of the trial wave function. 

We begin with the Hamiltonian H in its standard form f l2.1l )-( |?3D . Assume 
V{q) to be positive definite, and choose its minimum to be at g = 0, with 

V{q) > V{0) = 0. (^.1) 

Introduce a scale factor g'^ by writing 

Viq) = g\{q) {A.2) 

and correspondingly 

Thus, the Schroedinger equation (12.41) becomes 




where, as before, q denotes q\, q2, ■ ■ ■ , qN and V the corresponding gradient 
operator. Hence S{q) satisfies 

-^/(V^)' + ^^/V'^ + = E. (^.5) 
Considering the case of large g, we expand 

S{q) = So{q) + g-'Sr{q) + g-'S2{q) + • • • {A.6) 

and 

E = gEo + E^ + g-'E2 + • • • . (^.7) 
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Substituting (A. 6) - (A. 7) into (A. 5) and equating the coefficients of g ^ on 
both sides, we find 

V5o • V52 = ^ [V^^i - (V5i)2] - El, (A.8) 

etc. In this way, the second order partial differential equation (A. 5) is reduced 
to a series of first order partial differential equations (A. 8). The first of this 
set of equations can be written as 

l[VSo(«)f-«(5) = + . (A9) 

As noted in Ref.[15], this is precisely the Hamilton-Jacobi equation of a sin- 
gle particle with unit mass moving in a potential " — f (g)" in the N-dimensional 
g-space. Since g = is the maximum of the classical potential energy func- 
tion —v{q), for any point q ^ there is always a classical trajectory with a 
total energy 0+, which begins from q = and ends at the other point q ^ 0, 
with So{q) given by the corresponding classical action integral. Furthermore, 
So{q) increases along the direction of the trajectory, which can be extended 
beyond the selected point q 0, towards 00. At infinity, it is easy to see that 
<S'o(g) = 00, and therefore the corresponding wave amplitude e~^^°^'^'' is zero. 
To solve the second equation in (A. 8), we note that, in accordance with (A.l) 
- (A. 2) at g = 0, VSq oc v^{0) = 0. By requiring Si{q) to be analytic at 
g = 0, we determine 

Eo = ^(V'5o)at ,=0. (A.IO) 
It is convenient to consider the surface 

So{q) = constant; (^-H) 

its normal is along the corresponding classical trajectory passing through q. 
Characterize each classical trajectory by the ^o-value along the trajectory and 
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a set of N — 1 angular variables 



a = {ai{q),a2iq), • • • 



aN-i{q)) 



(A12) 



so that each a determines one classical trajectory with 



Vaj • VS'o = 0, 



(A13) 



where 



i = l,2,---,7V-l. (AM) 
(As an example, we note that as g — ^ 0, v{q) — ^ |Ecjfgf and therefore 

i 

So ^J2iOiqf. Consider the ellipsoidal surface 5*0 = constant. For 5*0 suf- 



ficiently small, each classical trajectory is normal to this ellipsoidal surface. 
A convenient choice of a could be simply any — 1 orthogonal parametric 
coordinates on the surface.) Each a designates one classical trajectory, and 
vice versa. Every {So, a) is mapped into a unique set (qi, q2, • • • , qN) with 
(So > by construction. In what follows, we regard the points in the (/-space as 
specified by the coordinates (So, a). Depending on the problem, the mapping 
(^1, q2, ■ ■ ■ , qN) (So, a) may or may not be one-to-one. We note that, for 
q near 0, different trajectories emanating from q = have to go along differ- 
ent directions, and therefore must associate with different a. Later on, as 
increases each different trajectory retains its initially different a-designation; 
consequently, using (So, a) as the primary coordinates, different trajectories 
never cross each other. The trouble-some complications of trajectory-crossing 
in ^--space is automatically resolved by using (So, a) as coordinates. Keeping a 
fixed, the set of first order partial differential equations can be further reduced 
to a set of first order ordinary differential equations, which are readily solvable, 
as we shall see. 



Write 



Si{q) = Si{So,a) 



the second line of (A. 8) becomes 




(A16) 
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and leads to, besides (A. 10), also 

S^iq) = SiiSo^a) = /-^[^V^^o - ^o], (^-17) 

where the integration is taken along the classical trajectory of constant a. 
Likewise, the third, fourth and other lines of (A. 8) lead to 

El = \[y^Si - (VSifU 5=0, (A.18) 
S2{1) = S2(So, a) = l ^^{-[V^Si - (VS^f] -E,}, (A.W) 

E2 = \[V^S2 - 2(VSl) • (VS2)]at ,=0, (^.20) 

5'o J c 1 

SM = S3(So, a) = / (^^{^[V'S^ - 2(VSi) • (VS2)] - B^}, (A.21) 

etc. These solutions give the convenient nornnalization convention at g = 0, 

5(0) = 

and 

e-^(^) = 1. (A.22) 

Remarks 

(i) As an example, consider an N-dimensional harmonic oscillator with 

V(q) = '^[ql + ql + --- + qly (A23) 

From (A. 2), one sees that the Hamilton-Jacobi equation (A. 9) is for a particle 
moving in a potential given by 

-v{q) = -\{ql + ql + ■ ■ ■ + ql). (A24) 
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Thus, for any point q the classical trajectory of interest is simply a straight 
line connecting the origin and the specific point, with the action 

So{q) = \{ql + ql + --- + q%). (A25) 
The corresponding energy is, in accordance with (A. 10), 

^0 = J. {A.26) 

By using (A. 8), one can readily show that Ei = E2 = ■ ■ ■ = and = 5*2 = 
• • • = 0. The result is the well known exact answer with the groundstate wave 
function for the Schroedinger equation (A. 4) given by 



and the corresponding energy 



E = ^g. {A.28) 



(ii) From this example, it is clear that the above expressions (A. 6) - (A. 8) 
are not the well-known WKB method. The new formalism uses —v{q) as the 
potential for the Hamilton-Jacobi equation, and its "classical" trajectory car- 
ries a 0+ energy; consequently, unlike the WKB method, there is no turning 
point along the classical trajectory, and the formalism is applicable to arbitrary 
dimensions. 

A. 2 Trial Function for the Quantum Double-well Potential 

To illustrate how to construct a trial function, consider the quartic potential 
( 11.31) in one dimension with degenerate minima. Set a = 1, we have 

V{x) = \g\x^ - If. (^.29) 
The corresponding Schroedinger equation is 



2dx^ 2 
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where, as before, ^(x) = e"^*^^^^ is the groundstate wave function and E 
its energy. Using the expansions (A. 6) - (A. 7) and following the steps (A. 8), 
(A. 10) and (A. 15) - (A. 21), we find the well-known perturbative series 

Sq{x) = ^{x-l)\x+2), Si{x)=\n^^, S2{x) = Yq~ 4^x'+\y ' "' 

{A.Sl) 

and 

Eo = l, Ei = ~ ^2 = "^' ^"^-^^^ 

Both expansions S = So+ g^^ Si+ g^"^ S2 + • • • and E = gEQ + Ei+g^^E2 + - ■ ■ 
are divergent, furthermore, at x = — 1 and for n > 1, each Sn{x) is infinite. 
The reflection x —x gives a corresponding asymptotic expansion Sn{x) — ^ 
Sn{—x), in which each x) is regular at x = —1, but singular at x = +1. 

As noted in section 1, for 5^ large, the first few terms of the pertur- 
bative series (with (A.Sl) for x positive and the corresponding expansion 
Sn{x) Sn{—x) for X negative) give a fairly good description of the true 
wave function i/j{x) whenever iIj{x) is large (i.e. for x near ±1). However, for 
X near zero, when ip{x) is exponentially small, the perturbative series becomes 
totally unreliable. This suggests the use of first few terms of the perturbative 
series for regions whenever tjj^x) is expected to be large. In regions where i/j^x) 
is exponentially small, simple interpolations by hand may already be adequate 
for a trial function, as we shall see. Since the quartic potential (A. 29) is even 
in X, so is the groundstate wave function; likewise, we require the trial function 
(f){x) also to satisfy (f){x) = (f){—x). At x = 0, we require 



(^).=o = 0. (A.33) 

To construct (f){x), we start with the first two functions So{x) and Si{x) in 
(A.Sl). Introduce, for x>0, 

1 I iZ/ 

and 

1 ~| 2C 
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In order to satisfy (A. 33), we define 

, . f + for 0<x<l 

^<^>^^<-^^ = {(l + ?^e'44.(i), for .^l 

Thus, by construction (A. 33) is satisfied. Furthermore, (\){x) is continuous 
everywhere, for x from — oo to oo, and so is its derivative. 
By differentiating 0+(a:) and we see that they satisfy 



and 
where 

and 

with, for X > 
where 



(r + y + w)0+ = #+ (A37) 



(r + y + /i)0 = #, (A38) 



1 



u{x) = , (^.39) 



h{x) = h{-x) (^.40) 
= u{x) + ^(x) (-4.41) 



v^ . 2ff (g-i)^^-^"^;^-;; , , forO<x<l ,,,,, 
for X > 1. 

Note that for g > \, g{x) is positive, and has a discontinuity at a: = 1. 
Furthermore, for x positive both u{x) and g{x) are decreasing functions of x. 
Therefore, hix) also satisfies ( 14.41) for a: = r > 0. 
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Appendix B. Proof of the Hierarchy Theorem 

In this Appendix, we give the proof of the Hierarchy Theorem stated in Sec- 
tion 4. The discussion follows very closely the one given for the one-dimensional 
case in Ref. 19. We first establish several lemmas applicable to both boundary 
conditions (]45D and (O): {A) /„(oo) = 1 and {B) /„(0) = 1. 
Lemma 1 For any pair fm{r) and fi{r) 

i) if at all r, ^ {^-^\ < then Sr,,+i > Si+i, (5.1) 

and 

ii) if at all r, f ] > then S^+i < Si+i. {B.2) 

Proof 



According to ( 14.91) 



^m+l[fm] = [h fm]- (-^.3) 

Also by definition 
Their difference gives 

(^m+l — ^;+l)[/m] = [{h — Sl+i) fm]- -5) 

From (B.3), 

0=[ih-8i+i) fi]. (5.6) 
Let n+i be defined by ( 14.211) : i.e., at r = n+i, 

hin+i) = Si+i. iB.7) 

Multiplying (B.5) by fi{ri+i) and (B.6) by /m(n+i) and taking their difference, 
we have 

fi{ri+i){Sm+i - Si+i)[fm] = [{h - Si+i) {fm{r)fi{ri+i) - fi{r)fm{ri+i))], 

iB.8) 
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in which the unsubscripted r acts as a dummy variable; thus [/m(^)] means 
[/„J and [/m(n+i)] means /^(n+i) • [1], etc. 

(i) If {fm{r)/fi{r)y < 0, then for r < n+i 

fm{r) fm{ri+i) (nQ\ 

Mr) /Kn+i) ■ ^ ^ ^ 

In addition, since in accordance with (14.41) and (B.7), h'{r) < and h{ri+i) = 
Si+i, we also have for r < n+i 

h{r) > Si+i. (5.10) 

Thus, the function inside the square bracket on the right hand side of (B.8) 
is positive for r < n+i. Also, the inequalities (B.9) and (B.IO) both reverse 
their signs for r > n+i. Consequently, the right hand side of (B.8) is positive 
definite, and so is its left side. Therefore, on account of (4.25A)-(4.25B), 
(B.l) holds. 

(ii) If {fm{r) / fi{r)y > 0, we see that for r < n+i, (B.9) reverses its sign 
but not (B.IO). A similar reversal of sign happens for r > n+i. Thus, the 
right hand side of (B.8) is now negative definite and therefore Sm+i < Si+i. 
Lemma 1 is proved. 

Let 

i = m {B.u) 

be a single valued differentiable function of ^ in the range between a and b 
with 

0<a<<^<6 (5.12) 

and with 

r/(a) > 0. (5.13) 

Lemma 2 . 

(i) The ratio 77/^ is a decreasing function of ^ for a < ^ < 6 if 

|<| at i = a iS.U) 

and ^ 

^ < for a<^<b. (5.15) 
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(ii) The ratio r]/^ is an increasing function of ^ for a < ^ < 6 if 

drj T] 



^>7 at ^ = a (5.16) 



and 



Proof Define 



^>0 for a<^<b. (5.17) 

(5.18) 



to be the Legendre transform L{^). We have 
and 

i (!) ^ 

Since (B.14) says that L{a) < and (B.15) says that ^ < for a < <^ < 6, 
these two conditions imply L(^) < for a < ^ < 6, which proves (i) in view 
of (B.20). The proof of (ii) is the same, but with inequalities reversed. 

Lemma 3 For any pair fm{r) and fi{r) 
(i) if over all r, 

d I fm{r)\ ^ d (Vm+l{r)\ ^ /noiN 

and (ii) if over all r, 

— —TTT- > then at all r, — — — ^ > 0. (5.22) 
Proof From and i KM . we have 



K+iir) = {h{r) - ^„,+i) x'{r) Ur) (5.23) 

and 

T^'iAr) = iKr) - x\r) Mr). (5.24) 
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Define 

^ = Vi+i{r) and r] = V^+,{r). (5.25) 

In any local region of r where Vi^^(r) ^ 0, we can regard t] = ?7(^) through 
Tiif) = r]{^{r)). Hence, we have 

*? ^ _ fs 26) 

. y -^"^' , (B.27) 



where 



and 



^ ^^i^'^ + f %V'| (5.28) 



fi V/. 



where 



(i) If {fm/ fiY < 0, from Lemma 1, we have 

^m+l > ^l+l- (5.30) 

From /i'(r) < and the definition of r^+i and n+i, given by (B.7), we have 

Tm+l < n+1, (5.31) 

h{rm+i) = and /i(n+i) = (5.32) 

We note that from ( l4l^ and ( |4^ ^^m+i(?^) and T>i+i{r) are both positive 



continuous functions of r, varying from at r = 0, 

P^+i(0) = P^+i(0) = (5.33) 

to at r = oo 

^^m+i(oo) = P^+i(oo) = (5.34) 
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with their maxima at r^+i for Vm+i{r) and n+i for Vi+i{r), since in accor- 
dance with (B.23)-(B.24) and (B.32), 

^^m+i(Wi) = and 2^i+i(n+i) = 0. (B.35) 

From (B.29)-(B.30), we see that R'{r) is always < 0. Furthermore, from 
(B.27), we also find that the function R{r) has a discontinuity at r = r/+i. At 
r = 0, -R(O) satisfies 

< R{0) = < 1. (5.36) 

h{0) - 8i+i 

As r increases from 0, R{r) decreases from R{0), through 

R{rm+i) = 0, (5.37) 

to — oo at r = r/+i— ; R{r) then switches to +oo at r = and continues 

to decrease as r increases from At r = oo, R{r) becomes 

i?(oo) = > 1. (B.38) 

It is convenient to divide the r-axis into three regions: 

(I) < r < rm+u 

(II) rm+i <r< n+i (5.39) 

(III) n+i < r. 

The signs of V^^^, and R' in these regions are summarized in Table 1. 

Assuming {fm/fiY < we shall show the validity of (B. 21), {Vm+i/Vi+i)' < 0, 
in each of these three regions. 



Table 1. The signs of V'^^^{r), P/+i(^)' H^) ~ ^m+i, h{r) - Si+i, 
R{r) and R'{r) in the three regions defined by (B.39), when 8m+i > £i+i- 



region 


^^Ui(0 






h{r) - Si+i 


Rir) 


R'ir) 


I 


> 


> 


> 


> 


> 


< 


II 


< 


> 


< 


> 


< 


< 


III 


< 


< 


< 


< 


> 


< 
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Since 

< hir) < S,^+i in II, (B.40) 

'^m+i{r) is decreasing and D/+i(r) is increasing; it is clear that (B.21) holds 
in II. 

In each of regions (I) and (III), we have R{r) > from (B.27) and R'{r) < 
from (B.29). Since {fm/ fi)' is always negative by the assumption in (B.21), 
both terms inside the big parenthesis of (B.28) are negative; hence the same 
(B.28) states that d?"ri/d^^ has the opposite sign from Vi^^. From the sign of 
V'lj^i listed in Table 1, we see that 

0<O in (I) (B.41) 

and 

^>0 in (III). (5.42) 

Within each region, 77 = ^^to+i(^) and ^ = P/+i(r) are both monotonic in r; 
therefore, 77 is a single-valued function of ^ and we can apply Lemma 2. In (I), 
at r = 0, both An+i(0) and 1^^+1(0) are according to (I4.20j) . but their ratio 
is given by 

^r..+ l(0) I^Ul(O) 



Therefore, 

(dri\ (t]\ 



(B.43) 



\dU \u ■ ^^-'^^ 

Furthermore, from (B.41), ^ < in (I), it follows from Lemma 2, case(i), 
the ratio r//^ is a decreasing function of ^. Since ^' = V'l^^ is > in (I) 
according to (B.24), we have 



d (T>m+i 



< in (I). (5.45) 



In (III), at r = 00, both 1)^+1(00) and 'D/+i(oo) are according to (B.34). 
Their ratio is 
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(5.46) 



which gives at r = oo, 



'dr]\ 



i 



(5.47) 



As r decreases from r = oo to r = n+i, from (B.42) we have ^ > in (III). 
It follows from Lemma 2, case (ii), rj/^ is an increasing function of ^. Since 

^' = V'l^i is < 0, because r > n+i, we have 



^(^]<0 in (III). 
dr V Vi+i J ^ ' 



{BA8) 



Thus, we prove case(i) of Lemma 3. Case(ii) of Lemma 3 follows from case 
(i) through the interchange of the subscripts m and /. Lemma 3 is then 
established. 

Lemma 4 Take any pair /m(^) and fi{r) 

(A) For the boundary condition /^(oo) = 1, if at all r, 



d_ I fmjr] 
dr \ fi{r) 

therefore, if at all r, 

_^ (fmir) 
dr \ fi{r] 



< then at all r, 



d ( fm+i{r] 



> then at all r, 



dr \ fi+iir] 



d ( fm+i{r] 



<0; 



^ A 4.U 4. II ^ I fm+lV) 

< then at all r, — ' 



dr V fi+i{r-] 



dr \ fi+i{r 

(B) For the boundary condition /,^(0) = 1, if at all r, 

_^ 1 fmjr] 
dr \ fi{r) 

therefore, if at all r, 

I fm{r) \ ^ Q ^^^^ ^ d_ ( fm+i{r] 

dr \ fi{r) ) ' dr\ fi+i{r) 

Proof Define 



> 0. 



>0; 



< 0. 



^ = fi+i{r) and 
From (14.151) we see that 

di fUAr) Vi^,{r] 



rj = fm+i{r). 



{BA9A) 



{B.50A) 



{BA9B) 



(B.bOB) 



(5.51) 



(5.52) 
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and 

d (df]\ Id (Vm+i{r] 



[A) In this case /n(oo) = 1 for all n. Thus, at r = oo, ^ = //+i(oo 
f} = /to+i(oo) = 1, and their ratio 



(5.53) 



= 1. (5.54) 



r=oo 



At the same point r = oo, in accordance with ( 14 .181 ). r>^+i(oo) = r>j„+i(oo) = 
0, but their ratio is, on account of h{oo) = and (B.l) of Lemma 1, 

i:>^+l(00) ^ ^^+i(00) ^ KOO) - 8^+1 ^ Srr^+l ^ ^ ,^ 

D;+i(oo) ^[+i(oo) h{oo) - 8i^i Si+i 

in which the last inequality follows from the same assumption, if {fm/ fiY < 0, 
shared by (B.l) of Lemma 1 and the present (B.49A) that we wish to prove. 
Thus, from (B.52), at r = oo 

'^^^ ^'"^'>(X\ =1. (B.56) 



r—oo 



As r decreases from oo to 0, ^ increases from /^+i(oo) = 1 to //+i(0) > 1, 
in accordance with (14.241) and (4.25A). On account of (B.21) of Lemma 3, 
we have {Vm+i/Vi+i)' < 0, which when combined with (B.53) and /^(r) < 
leads to 

d ( dfj 

diui, 



> 0. (5.57) 



Thus, by using (B.16)-(B.17) of Lemma 2, we have fj/i to be an increasing 
function of ^; i.e., 

d [fi 



Because 



> 0. (5.58) 



^(r^ 5 (I) ^"4(1,) 
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and f'l^i < 0, we find 

d ( fm+i\ d (t)^ 



< 0, (5.60) 



dr V fi+i J dr \^ 

which establishes (B.49A). Through the interchange of the subscripts m and 
/, we also obtain (B.50A). 

Next, we turn to Case (B) with the boundary condition /n(0) = 1 for all 
n. Therefore at r = 0, 

Furthermore from ( l4?T5l) and (4.16B), we also have /^+i(0) = //+i(0) = 
and Vm+i{0) = P/+i(0) = 0, with their ratio given by 

\dfi+Jr=o fUM -^i+m ^m(O) 
h{0) - 6i+i 

From (B.l) of Lemma 1, we see that if {fm/ fiY < 0, then S^+i > ^i+i and 

/ dfm+l\ 



therefore 



dfl+l J r=0 

I) 



< 1, (5.63) 
= 1. (5.64) 



Thus, 



dUr=Q \Ur=0 



Analogously to (B.18), define 



d^ Ji+i[r) 



therefore 



dL{r) _ ^,dX _ ^f^d_dv 
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From (B.21) of Lemma 3, we know that if {fm/ fi)' < then (D^+i/D^+i)' < 
0, which leads to 

dL{r] 



dr 

From (B.66), we have 



< 0. (5.68) 



and therefore at r = 0, because of (B.65), 

L(0) < 0. (5.70) 
Combining (B.68) and (B.70), we derive 

L{r) < for r > 0. (5.71) 
Multiplying (B.66) by //+i(r), we have 

/;+i(r)L(r) = /,+i(r)/,'„+i(r) - /,„+i(r)/;+i(r) 

Because //+i(r) and L{r) are both negative, it follows then 

which gives (B.49B) for Case (5), with the boundary condition /n(0) = 1. 
Interchanging the subscripts m and /, (B.49B) becomes (B.50B), and Lemma 
4 is established. 

We now turn to the proof of the theorem stated in (I4.26l )-( |434l) . 

Proof of the Hierarchy Theorem 
When n = 0, we have 

/o(r) = 1. (5.73) 
From (]4221)-(|424D, we find forn = 1 

/;(r)<0, (5.74) 
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and therefore 

{fi/foY < 0. (B.75) 
In Case (^4), by using (B.49A) and by setting m = 1 and / = 0, we derive 
(/2//1)' < 0; through induction, it follows then {fn+i/ fn)' < for all n. From 
Lemma 1, we also find Sn+i > Sn for all n. Thus, f l4.26l) -( j4~^^ are established. 

In Case {B), by using (B.75) and (B.49B), and setting m = 1 and / = 0, 
we find (/2//1)' > 0, which in turn leads to (/3//2)' < 0, and (]433)- 
i^Mj. Inequalities (lODD-JOS) now follow from (B.1)-(B.2) of Lemma 1. 
The Hierarchy Theorem is proved. 

Assuming that h{0) is finite, we have for any n 

0<Sn< h{0). 

Therefore, each of the monotonic sequences 

Si < 82 < S3 < ■ ■ ■ , in {A) 
Si < S3 < S^ < ■ ■ ■ , in (B) 

and 

S2> S4> Sq> ■ ■ ■ , in {B) 

converges to a finite limit S. By following the discussions in Section 5 of 
Ref.[18], one can show that each of the corresponding monotonic sequences 
of fn{f^) also converges to a finite limit f{r). The interchange of the limit 
n ^ 00 and the integrations in (4.16A) completes the proof that in Case (A) 
the limits S and /(r) satisfy 

/(r) = 1-2 fx-\y)dy j^x\z){Kz) - S)f{z)dz. {B.77A) 

As noted before, the convergence in Case (A) can hold for any large but finite 
and positive h{r), provided that h'{r) is negative for r > 0. In Case (B), 
a large h{r) may yield a negative /n(r), in violation of (4.25B) . Therefore, 
the convergence does depend on the smallness of h{r). One has to follow 
discussions similar to those given in Ref.[17] to ensure that the limits S and 
f{r) satisfy 

fir) = 1-2 l\~\y)dy l\\z){h{z) - S)f{z)dz. {B.llB) 



(B.76) 
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Appendix C. Comparison with Perturbative 

Expansion 

C.l A simple Model and Its Iterative Solutions 

Consider the one-dimensional problem in which V, U and h of (12 . II ) . (12.71 ) 
and (12.111) are given by 

oo \x\ > L + I 

V{x) = V{-x) = { for /<|x|<L + /, (C.l) 

-^g^ \x\ < /, 

cxD > L + /, 

U{x) = U{-x) = { for (C.2) 

^0 \x\< L + l 

and 

!0 \x\ > I, 

for (C.3) 

Let (f){x) and iIj{x) be the respective groundstates of 

W=-^0 + f/0 = ^o0 (C.4) 

and 

As a model, (C.5) is the Schroedinger equation that we wish to solve, and 
(C.4) is the equation that the trial function satisfies. The iterative equation 
(ng) becomes 

with i/jQ = Since U and V are both even in x, we need only to consider 

the region 

x>0. (C.7) 
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Set within x < L + I 

(j){x) = cospx, i.C.8) 

with 

2 

p{L + l) = ^ and Eo = ^. (C.9) 



Introduce, as in (12.251) . = We have in accordance with (l3.1j) -( |T7ll 



and therefore 



dD 

= (Jn{x) = {h{x) - Sn)(j)\x)fn-l{x), {C .11) 

ax 

oo 

J an{x)dx = (C.12) 



oo 9 

/ h(t)^fn-idx 

Sn = . (C13) 

/ (t)^fn-ldx 




One can verify that for n = 1 



g'- \ {l + '-^)[L + l-x-'-i!!^), l<x<L + l, 

and 



f(\ = ) fiiR) + ^[^-PiL + l-x)tanpx], l<x<L + l, 

~ [ /i(0) - - Si)xt8.npx, 0<x<l. ^ 

The constants fi{R) and /i(0) are related through the continuity of fi{x) at 
X = I. In order to derive the general solution for fn{x), it is useful to consider 
a limited region in x (either \x\ < I or \x\ > I) over which 

A„ = 2(/i(x) - = constant. (C.17) 
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From (C.lO)-(C.ll), we have 

~^(^°^^^^^) ^ (C.18) 

Write 

fn = Pn{x) + {i^npx)Qn{x). (C.19) 

We find, from (C.18), 

(f d 

—^Pn + 2p—Qn = -KPn-l 

dx^ dx 

and (C.20) 

Thus, P„ and Qn are both n*'^ order polynomials in x. Define 

P^^\{Pn^iQn) = {P;r . (0.21) 
From (C.20), P^ satisfies 

It is useful to introduce a set of n^^ order polynomials i?^, independent of A„. 
We require 

(£±2,,A)«.._«._, (C.23) 
with the boundary conditions (i) for n = 0, 

i?J(x) = 1 at ail X, (C.24) 

and (ii) for n > 0, at x = 

i?±(0) = 0. (0.25) 

The solution for Rn{x) can be further simplified by introducing a set of n^^ 
order polynomials i?„(a:|g) without the superscript + or — , but depending on 
a scaling parameter q, so that 



,x , 



Rnix\q) = q'^Rni-\l) (0.26) 
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where is an n*^ order polynomial of | with constant coefficients. The 

Rn{x\q) are related to Rn{x) defined by (C.23)-(C.25) through 



Rn{x\q) = Riix) when q = q+ = -\-^ 
Rnix\q) = R-{x) when q = q- = 



Writing 



we have 



correspondingly 



(0.27) 



u = - and Rni-\l) = rniu), (C.28) 
q q 

ro(w) = 1, ri{u) = u, r2{u) = \^ + w, 

Tziu) = lu^ + u^ + 2u, (0.29) 

D 

niu) = + + + 5w, etc. ; 



RQ[x\q) = 1, Riix\q) = qx, R2{x\q) = ^Q^x"^ + Q^x, 

Rs{x\q) = Iq^x^ + q'^x^ + 2g^x, (C.30) 
6 

115 
^4(a;|g) = t^Q'^x^ + x^^^^ + x/^^ + etc. 

In terms of Rn{x), the -P^(a:) that satisfies (C.22) is given by 

Pn{^) = ^AiA2---A„i?±(x) 

+c±A2A3---A,i?±_i(x) 
+c±A3---A,i?±_2(2^) + --- 

+c±_iAX(^) + 4 (C.31) 



where Ci , C2 , • • • are integration constants, with = (c^)* for all m. 
Thus for n = 1, since i?i(x|g) = qx, 

i^f (2:) = RM<t) = Rii^l ± ^) = (C.32) 
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Ai = 2{h{x) - Si) = 



According to (C.31), 

PHx) = Ia.RUx) +4 = dz^Ai + 4 iC.33) 
2 4p 

and from (C.21) 

Pi = p+ + Pf = c+ + cr 

and (C.34) 

= ii^p+ _ P-) = _|a^ + - cr). 

Since 

' -2£:i l<x<L + l, 

for (CSS) 
9^ -2ex 0<x<l, 

we have, in agreement with (C.16), Pi is a constant in / < a: < L + and a 
different constant in < x < /. Likewise, 

Qi = h constant \n I < x < L + I 

P 

and (0.36) 
1 

Qi = — ( Si)x + a different constant \n < x < I. 

These integration constants are determined by the continuity equations be- 
tween regions, leaving one final overall constant that is not determined. Like- 
wise, we can derive the solution fnix) for any n. 
If we impose the boundary condition at x = L -\- 1, 

UL + l) = l, (C.37) 

then in accordance with the Hierarchy Theorem 

8i<S2<--- (C.38) 

with 

]SmSn = S = E^-E (0.39) 

for any finite g^. 
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C.2 Analyticity of E{g^) 



The convergence of the perturbative series in g'^ depends on, in the connplex 
g"^ plane, the location of the nearest singularity g^ of E^g"^) to the origin. While 
the simple model discussed above can be readily solved, it is more complicated 
to locate the singularities of E{g'^). We first discuss the exact solution of (C.5). 
Again we need only consider the positive a:-axis. Note that the well-depth — ^ 
has a critical value — ^; for g'^ > g^, E < 0, and for g'^ < g^, E > 0. When 
g^ = gl ^ = and 



( L+l-x 
L 



1p =1pcOC < 



for 



cos QqX 



I < X < L + 1, 
< X <l 



with 



I 

T 



as required by the continuity of at x = I. When g^ < g^, 

' sm[k{L -\- 1 - x)] l<x<L-\-l, 

for 

< X <l 



COS qx 



with 
and 

When g^ > gl 

i/j oc 

with 



q tan ql = k cot kL, 
2 2^ ^ 



sinh[K(L — x)] 

cosqx 



for 



I < X < L -\- 1, 
< X <l 



q tan ql = k, coth k,L 



(0.40) 



(C.41) 



(0.42) 

(C.43) 
(0.44) 



(0.45) 



(0.46) 
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and 

S = -y = J(9^-fl^)<0. (CAT) 
In the limit L — ^ oo, go = 0, in accordance with (C.41), and therefore 

gl = 0. {CAS) 
For any however small, E is negative. Solution (C.45) becomes 



•0 oc < 

with 



e '^^ x> I, 

for (0.49) 
cosgx X < / 



and 

K, = qtanql. (C51) 

We examine the analyticity of E{g'^) for the case L = oo. the same analysis 
can be extended to arbitrary L, but it is more complex. The singularities of 
E = E{g'^) can be determined by setting = oo; i.e., ^ = 0, which leads 

"^^^ = 0. (0.52) 



to, on account of ^ = —Ik^, 



Because g'^ = + according to (C.50), (C.52) is equivalent to ^ + 1 = 0; 

^ = (0.53) 

Differentiating (C.51), we find 

— = tan ql + ql sec^ ql, (C. 54) 
dq 

which, when combined with (C.53), gives 

— — = tang/ + sec^ g/. (C55) 
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Since {q/n) = cot ql according to (C.51), (C.55) beconnes 

cot ql + tang/ + ql sec^ ql = 0; 

i.e., 

sec^ (//(cot ql + ql) = 0. 

When sec^ ql = 0, ql = ±ioo which is uninteresting for our purpose. Hence, 
for singularities at a finite g^, we require 

coi ql + ql = 0, (C.56) 

i.e., 

g/tang/ = -l. (CS?) 

Define 

z = g/. (C.58) 
Our problem is then reduced to the study of zeroes of 

ztanz+l = 0. (C.59) 

For each zero at z, there is a singularity of E{g^) in the connplex plane at 

g'^ = ti^ -\- q^ = q^ tan^ g/ + q^; i.e., at 

p^/^ = tan^ ^ + ^2 _ 1 + (^^_gQ^ 

(i) For z imaginary, we set z = ±iy. (C.59) becomes ytanhy = 1, which 
gives for the smallest \y\ = 1.2, z = ±1.2i and 

As will be shown below, this gives g'^. 

(ii) For z real, the nearest zeroes for (C.59) are z = ±2.8 and correspondingly 

g'^l'^ = 8.8 . (0.62) 

(iii) We shall now establish that in the complex g^-p\ane, the nearest singularity 
g"^ of -E(p2) to the origin is given by (C.61); i.e., 

gy = -0.44 . (C.63) 
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Consider the well known expansion 

c 

tan 2; = 2z 



00 



n=0 



For z near ^, we may write the righthand side as 



2z 

2z 
2z 



1 



+ E 

CO 

+ E 



1 



+ 



:n + i)47r4 



1 



1 4 „ °o 

- + (0--)+^ i: 



,-1 in + i)47r4 



+ 



(C.64) 



By using first the approximation 

tan z = 2z 



— (--—\ 



(C.65) 



and then the correction due to the successive remaining terms in (C.64), 
one can readily show that (C.61) gives the smallest singularity g1 in the 
^f^-plane. (Note that in this case, because L = 00, ^0 = in accordance 
with (C.48).) Thus, the perturbative series of E{g'^) is convergent only if 
(/2/2 < |(yf2|/2 ^ Q 44 On the other hand, the iterative series of Section C.l is 
convergent for any finite g^. 

C.3 A Comparative Study of Different Green's Functions 



(C.66) 



(i) Write the iterative equation (I2.19j) as 

(Hq - Eo)tljn{q) = Sn{q) 
where 

Sniq) = ihiq)-Sn)i^n-iiq). {cm) 

In one dimension, setting q = x and using (12 .251) and (l3.15n -( |TT6l ). the solu- 
tion il)n{x) is given by 



00 



(C.68) 
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where the Green's function Q is 

oo 

ix\g\z) = -20(x) / r\y)dyeiy - (0.69) 



with 




z> y, 

for (C.70) 



z<y, 



(x\g\z) = if z<x. (C.71) 

The function (xl^lz) is one of the Sturm-Liouville family of Green's func- 
tions[20]. On the other hand, in deriving the usual perturbative series one 
deals often with similar equations of the form 

{Ho - XHq) = v{q) iC.72) 

where A is a parameter. Let Un be the n^^ member of the ortho-normal set of 
real eigenfunctions of Hq with eigenvalue e„. We assume 

/ Un{q)um{q)d^q = 6nm = I ?' ^7^^: (C.73) 

Hence, 

{Ho - en)un{q) = (C.74) 
which reduces to (12.81) when n = 0; i.e., 

eo = Eq and uq = (f). (C.75) 

For A different from any of the eigenvalues Cn, the general solution of (C.72) 
is 

uiq)= Jiq\GiX)\qHq)d^q (C.76) 

where 

(<?|G(A)|g)^E ""^''^"'f^ . (C.77) 
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To simplify our notations, we assume all Un{q) to be real. When X = Eq = cq, 
(C.72) becomes 

iHo-eo)uiq) = viq), (C.78) 
same as (C.66). Since {Hq — eo)wo = 0, we have 

)v{q)d^q = 0. (C.79) 

The general solution of (C.78) is 

u{q) = j{q\G{\)\q)v{q)d''q + cu,{q) (C.80) 
with c an arbitrary constant and 

(q\G(e,)\q) ^ E . (C.81) 

In the sub-section (iii) below, we shall analyze the difference between Q and 
G, or G. Before that, we shall review the roles of G and G in the usual 
perturbative formulas.. 

(ii) In order to solve the Schroedinger equation (12.41) perturbatively, we begin 
with (EED; i.e., 

{H^-E^)i;{q) = {h-E)ilj{q). (C.82) 
Following the standard Brillouin-Wigner procedure, we require 

= / (t^{qmq)d''q = 1. (C.83) 



Define 

ij{q) = i^{q)-u,{q), {GM) 

which is orthogonal to UQ{q), because of (C.83). Eq.(C.82) leads to 

(Fo - ^o)V'(g) = {h- s){uo{q) + Hq)); (C.85) 

on account of (12.121) . this gives 

{Ho - E)ij{q) = {h- S)uo{q) + hi^{q). {GM) 

Multiply (C.86) on the left by Un and integrate. We have for n = 

S = J uoh{uo + i')d^q, (C.87) 
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the same as ( 12.161) : for n 7^ 

{en - E){n\ij) = J Unh{uo + ij)d''q (C.88) 

where 

(n|V^) = J UnH^q. {cm) 

Introduce the sequence 

^1, ^2, 4, ••• (C.90) 

with 

(0|^;) = for all / > 1. (C.91) 
When / = 1, we require for all n 7^ 

{en - E){n\iJi) = J Unhuod^'q. (C92) 

For / > 1, 

{en - E){n\i,i) = J Unh^i^id^'q. (C.93) 

Therefore, the sum 

i^ = Y.4^i {CM) 



satisfies (C.86). In compact form, 

i{q)= J{q\G{E)\q)S{q)d''q (C.95) 

where 

S{q) = hiP{q) + {h- 8)uo{q). {CM) 

The usual perturbative series can then be readily derived by first introducing 
a parameter e, defined by 

h{q) = eh{q). (C.97) 

The subsequent expansions of 

E = Eo + eE{l) + e'^E{2) + • • • (C.98) 

i^i = ei^i{l) + e%{2) + ■ ■ ■ (C.99) 
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and 

dGiEo) 



G{E) = G(Eo)+eE(l)- 



dEn 



lead to the usual perturbative series. 

(iii) According to (C.68)-(C.69), the Green's function Q satisfies 

(Ho - Eo){x\g\z) = 6{x - z), (C.lOl) 

whereas the function G{Eq) given by (C.81) satisfies, in a one-dinnensional 
problem, 

(^0 - E^){x\G[Eq)\z) = 5{x -z)- u^{x)u^{z) (0.102) 

with uq = (j) m accordance with (C.75). Thus Q and G are clearly different. 
As an explicit example, let 

= + (C.103) 

and 

x> R 

for < X < i?, (C.104) 

X < 0, 

similar to (C.2) and (C.4). The corresponding eigenfunction Un{x) and eigen- 
value Cn that satisfy (C.73)-(C.74) are 




Ur 



2 1 n + 1 

^ — sinknX and = -/c^, with /c^ = — ^tt, (C.105) 



andn = 0, 1, 2, For A different from any of the eigenvalues e„, G{\) 
of (C.77) in this one-dimensional problem is given by 

oo 1 

{X\G{\)\Z) = E vUn{x)Un{z) (ClOG) 
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where x and z can vary between and R. Introduce 

X+ = sinp(i? — x)^ X- = sinpx, (C.107) 

A = ip2 (C.108) 

and 

^ = (^)^+ " (^)^- = (^-109) 
The Green's function G{\) can also be written as 



(x\G(\)\y) = - 

Lu 



X-ix)x+iy) x<y, 

for (C.llO) 

x+{x)x-{y) x> y. 



One can readily verify that G{\) satisfies 

{x\G{\)\y) = {y\G{\)\x), (Clll) 

^~\&^ + U{^)){AG{\)\y) = 5{x - y) (ail2) 
and the boundary conditions 

(x|G(A)|0) = {x\G{\)\R) = 0. (ail3) 

As 

A-eo = i(^)^ (C.IM) 
the Green's function (^(eo) of (C.81) is given by 

^(eo) = lim (g(A) ^wo(y)wo(x)) 

A^eo V Co — A / 

= E — - — Un{y)un{x). (C115) 

Define 

= sin /co^c and 0(x) = cosA;o^ (C116) 
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with ko = ^. We find 

R . . . . 2 



{x\G\z) = -^(j){x)(f){z) {(f){x)z(f){z) + (f){z)x(l){x)] 



' for (ail7) 

4){z)(f){x) X > z. 
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The function ^ can also be expressed in terms of (f) and 0. We have 

= / Jf- (C.llS) 

and 

(^1^1^) = . I ^(^)^(^) - (f^iz)(i){x), x<z, (C.119) 
Thus, Q is related, but quite different from G and G. 
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